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The quantum Ising chain of length, L, which is separated into two parts by localized or ex- 
tended defects is considered at the critical point where scaling of the interface magnetization is 
non-universal. We measure the entanglement entropy between the two halves of the system in equi- 
librium, as well as after a quench, when the interaction at the interface is changed for time t > 0. 
For the localized defect the increase of the entropy with log L or with log t involves the same effective 
central charge, which is a continuous function of the strength of the defect. On the contrary for 
the extended defect the equilibrium entropy is saturated, but the non-equilibrium entropy has a 
logarithmic time-dependence the prefactor of which depends on the strength of the defect. 



I. INTRODUCTION 

Entanglement, quantum nonlocality and quantum cor- 
relations has became the subject of intensive research 
recently in different fields of physical: quantum informa- 
tion theory, condensed matter physics, quantum field the- 
ory, etc. For a quantum system which is divided into two 
parts, A and B, all information about entanglement is en- 
coded in the reduced density matrix: = Trg|^')(4'|, 
where |\E'} is a pure state of the complete system. The en- 
tanglement between A and B is conveniently measured by 
the von Neumann entropy Sa — — Tr^(/9^ logp^), which 
has been intensively studied in many-body systems, in 
particular in one dimension (Id). For a critical Id sys- 
tem (with periodic boundary conditions) the entropy is 
found to grow logarithmically with the length, L: 
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log L -I- ci 



(1) 



where L is the size of A or the size of the complete sys- 
tem, provided it is divided into two equal parts^>^>^. For 
conformally invariant systems the parameter in the pref- 
actor, c, is universal and given by the central charge of 
the conformal algebra. In the vicinity of the critical point 
where the correlation length is ^ <C i, the entropy is sat- 
urated and given by: 
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One is also interested in the time evolution of the 
entropy^ after changing the form of the interaction (quan- 
tum quench) at time i = 0. In the case of a local quench^ 
the interaction parameters arc modified in a restricted re- 
gion. For example measuring the entropy between A and 
B which are disconnected for i < but are joined to a 
closed chain with homogeneous couplings for i > at the 



critical point we observe a logarithmic increase in time. 
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If the complete system is open, i.e. there is one boundary 
point between A and B the prefactors in Eqs.((l][3]) are 
divided by a factor 2. 

Inhomogeneous interactions could modify the entan- 
glement properties of quantum spin chains. It has been 
shown that for rando m^i^°i^^i^^i^^i^'*i^^i^^ and aperiodioii 
couplings the prefactor in Eq. ([T]) is changed and involves 
the so called effective central charge, Cg//. On the other 
hand if the couplings vary linearly with the position an 
interface with a certain width is introduced, and in the 
expression of the entropy in Eq.© ^ is replaced by this 
length"'^^. 

If the inhomogeneities are centered at a few points 
("defects") they are not expected to modify the scal- 
ing form of the entropy, unless the defects are located 
at the interface. Indeed interface defects can modify the 
scaling form of the wavefunction in the vicinity of the 
junction^^, which in turn alter the entanglement entropy. 
The effect of a local interface defect. A, which measures 
the coupling between A and B has been investigated for 
XXZ and XX quantum spin chains2Si2ii^. For the an- 
tiferromagnetic XXZ chain the defect is a (marginally) 
relevant perturbation^'^, the defect renormalizes to a cut 
and the effective central charge approaches zerc^*'. On 
the contrary for the ferromagnetic XXZ chain the de- 
fect is a (marginally) irrelevant perturbation, the defect 
renormalizes to the homogeneous coupling and the effec- 
tive central charge approaches one. Finally, in the XX 
chain the defect is a marginal perturbation and the effec- 
tive central charge in Eq.([T|) is found^^ to depend on the 
strength of the defect, A. 

In the present paper we study the effect of interface de- 
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fects on the entanglement properties of critical quantum 
spin chains. Our approach differs from the previous ones 
in several respects. The system we consider is the quan- 
tum Ising chain and we study the problem with a local- 
ized as well as with an extended defect. The latter is real- 
ized by a smooth inhomogeneity in the couplings varying 
as ~ Ajx^ X being the distance from the interfaceS^^. 
Both perturbations are known to be marginal as far as 
the scaling behavior of the interface magnetization at the 
critical point is considered. We study the entropy both in 
equilibrium, as well as after a quench, when the interface 
couplings are modified for t > 0. The main goal of our 
investigations is to study possible relations i) between lo- 
cal critical scaling and the scaling of the entropy, and ii) 
between scaling form of equilibrium and non-equilibrium 
entropies like in Eqs. ([T]) and 

The structure of the paper is the following. The model, 
the type of defects, the interface critical behavior as well 
as the way of calculation of the entropy is described in 
Secini The localized and the extended defect problems 
are studied in Sec lIIII and llVi respectively. Our results 
are discussed in the final Section. Some technical details 
of the calculations are put in Appendices. 



II. MODELS AND METHOD 

We consider the quantum Ising chain defined by the 
Hamiltonian: 



Hf, = -"^ afa^^^ - /i V err 



1=1 



i=l 



(4) 



in terms of the Pauli-matrices, a^'' , at site i and with 
periodic boundary conditions, cr£^-^ — af. The quantum 
critical point of the system is given by^^ h = he = 1, 
where the bulk correlation function has a power-law de- 
cay for large L: 



G{L)^ 
with r/ = 7^0 = 1/4. 
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which is a continuous function of the strength of the de- 
fect, A. Note that in the special cases, A = and A = 1 
we recover the known results for surface and bulk corre- 
lations, respectively. 

More generally, one can consider the defect at posi- 
tion L of different strength, say A'. For A' = 1 and for 
A' = there is one defect in the closed or open chain, 
respectively. The defect exponent is then modified to 
{■qd + l/4)/2 and {-qd + l)/2, respectivelj*i^. 



B. Extended defect 

The extended defect is defined by a smooth 
inhomogeneity^^: 
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and this perturbation is put symmetrically in the two 
parts of the lattice. In Eq.® we have used a shift, 5 — 
0(1), in order to avoid singularities. 

The local critical behavior of this system is different 
for A < and for v4 > 0. For weakened local couplings, 
A < 0, correlations between two defect spins, Gd{L), has 
an algebraic decay with an A dependent defect exponent: 
Vext{A) — 1 ~ A. At the same time correlations between 
two bulk spins, G{L), involves an exponent rj which also 
depends on A: rj — l/A — A. This unexpected varia- 
tion of the exponent with A is due to the fact that the 
interface coupling at i = (as well as at i = L/2) renor- 
malizes to zero as Jq{L) ~ L^. This fact explains also 
the observation that the decay exponent of the end-to- 
end correlations for two decoupled chains is just rjextiA), 
at least for ^ < 0. 

For enhanced local couplings, A > 0, the defect-defect 
correlations, Gd{L), approach a finite limiting value, m^, 
so that the interface stays ordered at the bulk critical 
point. The connected correlation function, Gd{L) — m^, 
decays to zero algebraically with an exponent ry^^j (A) = 
2 A. In this case correlations between bulk spins, G{L), 
involves the pure exponent, 770 = 1/4. 



A. Localized defect 

A localized defect is defined by the perturbation: 

Vioc = (1 - A)(a2/2a2/2+i + <^l<^l+i) (6) 

so that the complete Hamiltonian is given by Hq + Vioc- 
This perturbation does not modify the decay of the bulk 
correlations in Eq. (O , however the interface or defect cor- 
relations: 

Gd{L) = {Q\a^,al^,\Q) ^ L-^- , (7) 
involve a new exponen t^^i^^ : 

?7d = ?7,oc(A) = \ arctan2(l/A) (8) 



C. Calculation of the entropy 

Calculation of the entropy of the quantum Ising chain 
in the equilibrium case is described in detail in several 
papers^ii^i^. The nonequilibrium entropy in the ho- 
mogeneous case and in the thermodynamic limit is cal- 
culated in^, whereas for inhomogeneous couplings the 
method is described in the Appendix [A] Here we briefiy 
recapitulate the main steps of the calculation and de- 
scribe the technical steps needed in the numerical calcu- 
lation. 

The first step is to transform the Hamiltonian of the 
quantum Ising chain in terms of free fermions^^. Numeri- 
cally, this step necessitates the diagonalization of an L x L 
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symmetric matrix. In the second step we calculate the 
reduced density matrix, which can be reconstructed from 
the correlation matrix in the free fermionic basis^. The 
entanglement entropy is calculated then from the eigen- 
values of the reduced correlation matrix., If A has £ sites 
(in our case we use £ = L/2) this second step requires the 
diagonalization of a, £ x £ symmetric matrix, if one works 
with Dirac fermion&i^ or a 2^ x 2£ skew-symmetric matrix 
if the calculation is performed with Majorana fermions^. 
At this step we have calculated the so called single-copy 
entanglemenli^, too, which is defined through the largest 
eigenvalue of the density matrix, wi, as Si — — log wi . 
Si is obtained also from the eigenvalues of the correlation 
matrix. For homogeneous chains 5*1 (L) is known to have 
a logarithmic size-dependence: 
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with a prefactor: k — c/2. 

In the third step of the calculation we consider the 
time evolution of the entropy after a quench. For this 
we should calculate the time evolution of the correlation 
matrij*^, which can be done in the Majorana fermion 
basis. Each matrix-elements is obtained through op- 
erations, which will result in the increase of the compu- 
tational time accordingly. 

If the entropy has a logarithmic dependence, either as 
a function of the size (see Eq.(IT])) or the time (see Eq.®), 
then we have calculated the prefactors through two-point 
fit and in this way effective central charges are obtained 
as: 



(L) = 3[5(2i)-5(i)]/log2 



(11) 



Similarly in the nonequilibrium case we calculated the 
prefactors as: 

p{L) = [5(2i,<-i/2)-5(i,< = i/4)]/log2. (12) 

From this series of results an estimate of Ce/ / or p is ob- 
tained by sequence extrapolation methods, such as by the 
Bulirsch-Stoer algorithm^^. In the numerical calculation 
we used finite systems up to L = 1024 for the equilib- 
rium entropy and up to L = 512 for the non-equilibrium 
entropy. 



III. CHAIN WITH LOCALIZED DEFECTS 

Having two symmetrically placed defects of strength A 
in the quantum Ising chain, as given by the Hamiltonian 
Ho + Vext in Eqs.Q and Eq.® we have calculated the 
entanglement entropy between two halves of the system 
for different lengths, L. The results are shown in Fig[T] 
as a function of logL. 

For large L the curves approach straight lines with dif- 
ferent A-dependent slopes. We have calculated effective 
central charges by two point fits, see Eg. pTjl . and their 
extrapolated values are put in Fig. [21 These are contin- 



FIG. 1: (Color online) Entanglement entropy of the quantum 
Ising chain with different strength of the defect as a function 
of the logarithm of the length. Inset: the same for the single- 
copy entanglement. 
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FIG. 2: (Color online) Effective central charge, Ce//(A) (full 
line - green), and the parameter, Ke//(A) (broken line - red), 
of the single-copy entanglement of the quantum Ising chain 
as a function of the strength of the localized defect, A. For 
small A the leading behaviors obtained by perturbational cal- 
culation are indicated by dotted lines. 



uous function of the strength of the defect and vary from 
to 1/2 as A tuned from to 1. The numerical data 
are consistent with the relation, Ce//(A) = Ce//(1/A), 
which symmetry holds also for the local magnetization 
exponent in Eq.®. For small A we have calculated S{L) 
perturbatively, the calculation is presented in Appendix 
IB] We have obtained that in leading order of A^, S{L) 
has a logarithmic L-dependence and the effective central 
charge is given by: 

Ce//(A) = 6A' (^-^{1 - In A^) + b^+ OiA^) . (13) 

Here b = 0.062180(2) is a numerically calculated con- 
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stant. The perturbative result in Eq. (fT3|) is also shown in 
Figl2] together with the numerical data. It is interesting 
to note that in Eq. (fT3| there is a logarithmic correction 
to the leading behavior. 

We have also calculated the single-copy entanglement, 
Si, which is found to be in the form of Eq. pH)) . however 
with A dependent prefactors, Kg// (A), the extrapolated 
values of which are plotted in FiglH too. In the range 
< A < 1, Keff{^) is seen to vary between and 1/4 
and in the small A limit we have (see Appendix IB|) : 

«e//(A) = — + 0(A4). (14) 

TT 

The ratio Ke//(A)/ce//(A) varies between and 1/2, 
thus the conformal result, n/c— 1/2 is valid only in the 
homogeneous system. 

Our results can be generalized if there is one defect 
in the system. For the closed chain (with A' = 1) the 
effective central charge is changed to (ce//(A) + l/2)/2, 
whereas for an open chain (with A' — 0) the prefactor of 
the logarithm is just the half as for two defects. 
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FIG. 3: (Color online) Time evolution of the entropy after 
a quench from two disconnected chains. The final system is 
closed and homogeneous (a) or contains one (b) or two (c) 
defects of strength A = 1/2. The quench to a homogeneous 
open chain is given by (d). The analytical formulae for the 
homogeneous closed Ea. (|15p and open Ea. (|16p chains are also 
shown by full lines. 



A. Time evolution after the quench 

Here we calculate the time evolution of the entropy by 
starting with two disconnected parts, i.e. with A = 
(and A' = 0) for i < 0, and connecting them by one 
or two defects with A > for i > 0. In a finite chain of 
length L the entropy has a periodic time-dependence, the 
period of which is L/2, if the final chain is closed (A' = A 
or A' = 1.) and the period is L, if the final chain is open 
(A' = 0.). This is illustrated in Figl3] for a chain with 
L — 128 sites and with different type of defects. 

If we perform the quench to the homogeneous chain, 
the numerical data are very well fitted by the formula: 
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(16) 



for an open chain with c — 1/2, which are also shown in 
FigEl 

In the general case, when the closed chain for i > 
contain one or two defects of strengths < A < 1, the 
numerical results in Fig[3]are compatible with the scaling 
form (for closed chains): 



est 



(17) 



where the scaling function, f'^^{y) > 0, is periodic with 
period 1/2 and for small argument it behaves as f'^\y) ~ 
y. For an open chain with one defect the prefactor in 



Eq. p?|) is changed to its half and the scaling function, 
f°^{y) > 0, is periodic with period 1. Then for t ^ L we 
expect the asymptotic behavior: 



5°P(t) ^ logi + cst 



(18) 



which is a generalization of Eq.Q. This relation is 




FIG. 4: (Color online) Time evolution of the entropy after a 
quench from two disconnected chains into an open chain with 
one defect having different strengths A = 1. (homogeneous), 
A = .75, A = .5 and A = .25, from the top to the buttom. In 
the logt scale the initial part of the curves for L = 128 and 
L — 256 are close to the indicated straight lines having the 
slope Ce//(A)/3, which is calculated from the scaling of the 
equilibrium entropy (see Fig[2J. 

checked in FigH] in which the time evolution of the en- 
tropy is shown as a function of log t for different values 
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of A. Indeed the starting part of the curves are well 
described by straight lines the slope of which is compat- 
ible with Ce//(A)/3, as calculated from the equilibrium 
entropy and given in Fig[21 

We have repeated the calculation by considering an- 
other type of quench: for t < the system contains a 
pair of defects of strength < A < 1, which is changed 
to the homogeneous coupling, i.e. A = 1 for < > 0. In this 
case the scaling form in Eg. p?)) is still applicable, how- 
ever with a different effective central charge, Cg^j(A), 
which is a decreasing function of A. We have checked 
that Ce//(A) -I- Cgyy(A) = c(A) > 1/2, for example 
c(.25) = 0.527(3), c(.5) = 0.545(5) and 5(.75) = 0.520(5). 



IV. CHAIN WITH EXTENDED DEFECTS 

We have calculated the entanglement entropy between 
two halves of the quantum Ising chain which contains a 
pair of extended defect, as described by the Hamiltonian 
Ha + Hext in Eqs.lg]) and For different lengths of the 
chain, L, the entropy is plotted in Fig|5]as a function of 
the strength of the defect, A. 
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FIG. 5: (Color online) The entanglement entropy of the quan- 
tum Ising chain with two symmetrically placed extended de- 
fects as a function of the strength of the defect, A, for differ- 
ent finite systems. The conjectured asymptotic behavior in 
Ea. (l20|l is indicated by the full red curve. Inset: Scaling plot 
of the entanglement entropy around its maxima {exp{S — Sm)) 
in terms of the combination {A — A,n){log2 L — 3), see the text. 



For a given L the entropy has a maximum close to 
A = 0, which corresponds to the critical pure system, 
whereas for large negative (positive) A the entropy ap- 
proaches the limiting value (log 2), which is the same in 
the fully disordered (ordered) phase of the pure system. 
For intermediate values oi A ^ the entropy with in- 
creasing L seems to be saturated. This follows from the 
observation, that the effective central charges obtained 
through Eq. Jlip approach zero even for a small value of 
|^|. 



The finite-size dependence of the entropy can be ex- 
plained if we take into account the critical scaling be- 
havior at the interface, which is valid for large enough 
lengths, L > L. As described in Sec|TTl for A < the de- 
fect renormalizes to a cut, whereas for A > it becomes 
ordered, which is in agreement with the behavior of the 
entropy in Fig[51 The microscopic length-scale, L, can 
be estimated in the ^ < regime from the value of the 
renormalized connecting coupling: Jo(-^)Il=l ^ L"'"*', 
which should be in the order of 0(1), say e~^, with 7 > 0. 
From this we obtain for the microscopic length: 



exp 



(1^1 



(19) 



A similar expression can be derived in the regime ^ > 0, 
too. Note that L has a very fast increase with decreas- 
ing \A\ and it is divergent in the homogeneous system. 
The microscopic length, L, sets in a length scale for the 
entanglement, too, and in the limit L <^ L the entangle- 
ment entropy is obtained from Eq. ([2|) by replacing ^ with 
L. so that 



S{A) c^- log L 



- est ~ - 
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We have tried to fit the extrapolated numerical values 
of S{A) by this formula in Figl5l which is found to be 
reasonable with 7 ~ 0.87. 

We have also analyzed the finite-size scaling behavior 
of the entropy close to its maximum, which is located at 
A„i — AmiL) and has a value Sm = Sm{L). In order 
to shift the maximum to the same position for all L we 
have considered the difference, AS{L) = S{L) — Sm as 
a function of A A = A — Am- In terms of the scaling 
variable AA(logL— 3) the curves for different L are scaled 
together, as illustrated in the inset of FigO Here the 
scaling collapse is very good for A < 0, whereas for A > 
the somewhat less perfect collapse is probably due to the 
presence of interface order in the system. Since Am{L) — > 
for large L the scaling combination used in the inset 
of FiglS] is compatible with the microscopic length-scale 
defined in Eq.dll 



A. Time evolution after the quench 

We have measured the time evolution of the entropy, 
if the system contains a pair of extended defects for 
t < 0, which is removed for t > and we are left 
with the homogeneous closed chain. As shown in FiglS] 
the entropy for t <C L has a logarithmic increase in 
time and the prefactor, p{A), which seems to have the 
symmetry: p{A) = p(—A) is increasing with \A\. We 
have calculated two-point fits for the prefactors, see in 
Ea. (fl2l) . the extrapolated values of which are plotted 
in the inset of FiglG] The prefactor has its minimum 
around ^ = 0, which is close to 2c/3 = 1/3, whereas 
for large \A\ we have approximately: p{A) « A/2. The 
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FIG. 6: (Color online) Time evolution of the entropy after 
removing a pair of extended defect of various strength, A, 
from the critical quantum Ising chain (for t > 8. from the 
top to the bottom A=-8.;-6.;-4.;4.;-2.;2.). In a logt plot the 
initial part of the curves are described by straight lines with 
A dependent slopes. The extrapolated values of the slopes 
as a function of A are given in the inset together with an 
interpolation curve, see the text. 



measured points can be well interpolated by the curve: 
p{A) = [1/9 + \A\/6 + which is also indicated 

in the inset of FiglHl 

To explain the observed behavior of the entropy we 
have to take into account two different effects of the ex- 
tended defect. First, from local scaling consideration in 
Sec|ll]the defect renormalizes to a cut, thus - for small 
\A\ - the quench is made from two disconnected parts to 
a closed chain, and according to Eq.([T5]) the prefactor 
is 2c/3 — 1/3 in agreement with the measured limiting 
values. The second effect of the extended defect is to pro- 
duce an increase of the entropy through inhomogeneous 
quench, since the interaction in the chain at position x 
differs from the critical value by J{x) — 1. According to 
the argument in Refi^ pairs of quasiparticles are emitted 
at different points of the chain and they will contribute 
to the entanglement at later time, when reach the two 
parts of the system. The density of quasiparticles, a{x) 
depends on the distance from the critical point at the 
given position'^*, and for a small perturbation it is given 
hy^ a ~ \J{x) — l|/2. The increase of the entropy is 
obtained by integrating over the contributions: 

S{t)-SiO) = ^j'^dxa[Jix)], (21) 

which for the smooth inhomogeneity in Eq.(I5]) and for 
t <^ L results in: S{t) — 5(0) « A/2\ogt, in agreement 
with the numerical results. 



V. DISCUSSION 

In this paper we have considered the quantum Ising 
chain with critical couplings which is separated into two 
parts by localized or extended defects. In both cases the 
scaling behavior of the interface magnetization is non- 
universal (the scaling exponents depend on the strength 
of the defect) and we asked the question, how this fact 
is reflected in the entanglement properties of the sys- 
tem. For the localized defect both the equilibrium and 
the nonequilibrium entropy is found to be characterized 
by the same effective central charge, the value of which 
depends on the strength of the defect. In this case scaling 
of the equilibrium and the nonequilibrium entropy can be 
cast into the same form, see Ea. ([T7|) . 

The situation is found completely different for the ex- 
tended defect, in which case the equilibrium entropy is 
saturated, although the magnetization correlations (also 
at the interface) are long ranged and thus the correspond- 
ing correlation length is divergent. In this case the entan- 
glement is related to another, finite microscopic length, 
as given in Eq. (fT9)l . This is due to the fact that the two 
parts of the system at the interface are asymptotically 
separated and the wave-function become localized. We 
can thus conclude that for an extended defect the corre- 
lation length and the entanglement length have different 
scaling properties. As far as the nonequilibrium entropy 
of this system is concerned it is shown to have a logarith- 
mic ^-dependence, the prefactor of which is the result of 
two effects: the asymptotic cut and the inhomogeneous 
quench. 

Our results for the localized defect are related to sim- 
ilar studies for the XX-cYv&m^^-^^ . Since the entropy 
of the quantum Ising chain and that of the XX-cYiaJm 
are exactly related^ the same is true for the central 
charges, too. For example from Ce//(A) in Fig[2] we 
obtain the effective central charge in the closed XX- 
chain with one interface defect of strength < = A as 
c^/j{t) = Ce//(A)-hl/2, as given in Fig. 5 of Ref,^. Here 
we comment on the observation in Refi^ that the small 
A dependence of the effective central charge is given by: 
Ce//(A) ^ , with an exponent, 5 « 1.8. According to 
our perturbative calculation in Eq. (jl3[) the true exponent 
is 6 — 2, however with a multiplicative logarithmic cor- 
rection term-^. Similar behavior is expected to hold in 
higher dimensional gapless fermionic systems with weak 
links^. The result in Eq.((T3]) can be compared with a 
bosonization study of the continuum version of the XX- 
chain^, in which the impurity (defect) contribution to 
the entanglement entropy is found to scale as: 

dS=^yhHn^ (^"^""t) ' ^^^^ 

where y is the strength of the impurity potential and e 
is a small-distance cutoff. The leading ln^(L/e) term in 
Eg. ip^ is not compatible with our lattice result. 

As far as marginal extended defects are concerned our 
results for the nonequilibrium entropy are expected to 
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be generic for another critical quantum spin chains, too. 
From scahng theory it is known that an extended de- 
fect in the form, (A/x)", x being the distance from the 
center of the defect, is a marginal perturbation, provided 
uj = l/v, V being the correlation length exponent^^. Now 
estimating the nonequilibrium entropy we repeat the ar- 
gument at the end of Sec lIVi where the density of quasi- 
particles is expected to scale with the value of the local 
gap2^: a{x) ^ \J{x) — l]" ^ A/x. Integrating the con- 
tributions in time, see Eq.([2T|). for large A leads to the 
behavior, S{t) — <S(0) ~ Alogt, as for the quantum Ising 
chain. 
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tional Research Fund under grant No OTKA T048721, 
K62588, K75324 and M045596 and by a German- 
Hungarian exchange program (DFG-MTA). We are 
grateful to I. Peschel and H. Rieger for useful discussions. 



APPENDIX A: TIME EVOLUTION OF THE 
ENTROPY FOR QUADRATIC FERMIONIC 
SYSTEMS 

Let us consider a general Hamiltonian, Ti, which is 
quadratic in terms of fermion creation, c^, and annihila- 
tion, Ck, operators and which is given for i > as: 



k,l = l 



clAkici + i (cis, 



klCi 



(Al) 



Here Aki = {A)ki = Aik and Bui = (B)^; = -Bik are real 
numbers and k, I are the sites of a lattice. In the initial 
state, i.e. for f < the parameters of the Hamiltonian 



are different, say, A^^i and B^;'' and the ground state of 
the initial Hamiltonian, t6^\ is denoted by |V'o}- The 
system is divided into two parts: A consists of points 
k ~ 1,2, ... ,£ and B of the rest of the system. 

For one dimensional spin models, such as the quantum 
Ising chain, the Pauli spin-operators, af'^'^, are related 
to the fermionic operators as^: 

[](-a,)Vf = cj + Q = 2t, = {~iy-^a2i-i 
j<i 

l[{-<J,yiaf = cj -ci = ^i = {-iy-^a2i . (A2) 

3<l 

Here, in the last two equations, at a given site, I, we 
define two Clifford operators, 2t; and 03/, as well as two 
Majorana fermion operators, 0,21-1 and 0,21. The commu- 
tation relations for these new set of operators are: 

21? = 1, 05? = -1, = -*Bfc2l/, 

2li2lfc = -2tfc2ti, Q3i<Bfc = -OSfeS,, //fc (A3) 

for the Clifford operators and 



(0) 



for the Majorana fermion operators. 



The first step of the calculation is to diagonalize both 
, H^") and H, which can be made by the same type of 
canonical transformation. For simplicity we work here 
with Ti, for Ti^*^^ the analogous results are denoted by 
a superscript '^^^ . The new set of fermion operators are 
given by the combination: 



Vk 



= '^{gkici + hkicj) 



(A5) 



where gki and hki are real numbers and the Hamiltonian 
assumes the diagonal form: 



n = J2Akvl 



VkVk 



const. 



(A6) 



fc=i 



Here the energy of the free- fermionic modes, A^,, are given 
by the solution of the eigenvalue equations: 



(A-B)(A + B)*fe = Al^, 



^l^k 



(A + B)(A-B)M>fe 
and the components of the eigenvectors are: 



*fe(i) = gki 
*fe(«) = gkt - 



^ hki 
hki ■ 



(A7) 



(AS) 



In the second .step we consider the subsystem A and cal- 
culate its reduced density matrix, p^, which can be re- 
constructed from the time-dependent reduced correlation 
matrix of the Majorana operators: 



{V'o|am(i)a„(i)|'0o) = S, 



m,n "I" ^(r^ ) 



(A9) 



m,n = 1,2, . . . ,2£. Here {Tf)mn is a skew-symmetric (or 
antisymmetric) matrix which is transformed by an or- 
thogonal transformation, Q, into a block-diagonal form: 



■ i/i 

-1^1 



-^2 



QPfQ^ 



. l^r 

- Vr 



(AlO) 



thus the eigenvalues of Vf are ±iVr, r = 1,2, . . . In 
this representation the reduced density matrix is the di- 
rect product of £ uncorrelated modes: pg — ^^.^iPr, 
where pr has eigenvalues (1 ± Vr)l2. The entanglement 
entropy is then given by the sum of binary entropies: 



(A4) ^l(^) 



/ 1 + f ,. , 1 + f 1- Vr . 1 - V, 

— — log2 — \ — log2 



2 

(All) 
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In the third step we calculate the time-dependent corre- 
lation matrix and work in terms of the Clifford-operators 
the time evolution of which are given bj^: 

^i{t) = J2 [(^i^fe)t2lfc + {AiBkh'Bk] 

k 

= [{BiAk)t'^k + {BiBk)t'&k\ (A12) 

k 

Here the time-dependent contractions are: 

{AiAk)t = ^cos(Agi)$,(0$g(fc) , 

q 

{AiBk)t = (BfeAz)t =z^sin(A,t)$,(0*g(fc) , 
{BiBk)t = ^cos(A,t)*,(0*g(fc) . (A13) 

Note, that in Eq. (jAl3P the free-fermionic quantities 
are related to the Hamiltonian 7i, which governs the 
time evolution in the system for t > 0. The matrix- 
elements of the time-dependent Clifford operators, such 



APPENDIX B: PERTURBATIVE CALCULATION 
OF THE ENTROPY FOR LOCALIZED DEFECTS 

Here we consider the homogeneous quantum Ising 
chain with one single defect coupling of strength, A, 
between spins at £ and £ + 1, but the two ends of the 
chain at 1 and L are free. The unperturbed system (with 
A = 0) consists of two separated chains: A with sites 
i = 1, 2, ... £ and B with sites i ^ £ + 1,£ + 2, . . . L, and 
the length of B is denoted hy £' — L — £. The reduced 
density matrix of „4, denoted by pe, is calculated pertur- 
batively in Ref.-^^-. The matrix-elements in leading order, 
i.e. up to O(A^) are given by: 

B 

{ipf\p,\cpf) = pS,j) ^ zY,c{t,k)c*ij,k) , (Bl) 

k 



as {ipo\^i{t)^k{t)\iljQ) , involve the following ground-state 
expectation values: 

{i,o\^k,^kM = (Vo|«fc,a,jVo) = Gi°iAi4) 

Here the first equations follow from the commutation 
rules in Eq. (|A3p . whereas the static correlation matrix 

Gil = -E*r(fciX^(M, (A15) 

which is calculated with the initial Hamiltonian, 7^*^°-'. 
Now one can go on and calculate the time-dependent 
correlation matrix of the Clifford operators, which is then 
transformed by Eq (|A2p for Majorana operators so that 
finally one obtains the matrix-elements of in Eq. (|A9p 
as: 



(A16) 

I 

with c(0,0) = 1 and 

ci^,k) = ■ (B2) 

Here |(y9^) (|</jf )) is the z-the (fc-th) eigenstate of the un- 
perturbed system A (B) with excitation energy: (E^) 
and Z is a normalization constant, so that pe{i, i) — 
1. In Ea. (|B2[) the matrix-elements of the surface mag- 
netization operators af of system A and crf_|_i of system 
B are non-zero, if the excited states, and ) in 

the fermionic representation contain just one fermion. 
Then using the notation of AppO the non-vanishing 
matrix-elements are given by: {Vi'\'^f\'-P'o) = ^i'i^)^ 
i = with Et = A^ and laf+il^p^) = vE-f (1), 

k — . . .£' with E^ = A^. Using the exact solution of 



- 2/-l,2m-l 



^ 2i-l,2m 



2l.2m-l 



2i,2n 



J2 G^^X^{AiBk,)t{A„Mt-Y. Gil(^;^'c,)*(A„i?fe,), 
^ G^^^^{AA,,),{B,Mt~Y. G'^^,.^{AiB^,)t{B„Mt 

ki ,A,'2 ki ,k2 

GfXii^-'n^ki)t{BiBk2)t - E G^,^^,^^{AmBk^)t{BiAk:,)t 

k\,k2 ki,k2 

J2 GfX,{BiAkMBrnBk,)t - J2 G'^kL{BiBu,)t{B„Mt 



fei ,k2 



ki ,k-2 
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the open quantum Ising chain with homogeneous critical 
couplings of length i (see Eqs.(A.5) and (A. 6) of Ref.^^) 
we obtain in leading order: 



E 



cos Ui COS aj 
(2^+l)(2r + l) 

cos^ 



^ (sinofj + sin/3fc)(sinaj + sin/3fe) 



(B3) 



1,2, 



, £, with Of,- 



2i-l J O _ TT 2k-l 

2 2£+l /^fc 2 2£' + l • 



In 



the following we analyse the consequences of this expres- 
sion in the hmits: f > 1 and £'/£ > 1. 

The leading eigenvalue of the reduced density matrix 



u;i = Z = 1 - (ai + alogi) , 



(B4) 



where the prefactor of the logarithm is a = tt^^, and oi 
is a constant of 0(1). From Eq.(jB4| we obtain for the 
single- copy entanglement; 



Si 



A^(-log^ 



est. 



(B5) 



which is proportional to log£. For two symmetrically 
placed defects, as studied numerically in Sec lIIII the pref- 
actor is 2A^/7r^, from which we obtain the value of 
Kg// (A) as given in Eg. pi)) . 

To obtain the entanglement entropy in leading order 
one should calculate the other i = 2,3, . . . £ + 1 eigenval- 
ues of the reduced density matrix, which are expressed 
as: Wi — A^Ei, so that the entropy is given by: 



5(A) = -?«ilogu;i-A2^e,(loge,-KlogA2) + 0(A4) . 

i=2 

(B6) 

Here the correction term is evaluated numerically and 
we have obtained J2^=2 €i — bi + b log £ and the 
prefactor of the logarithm is b = 0.062180(2). Putting 
this and wi from Eq. (jB4p into Eq. (|B6p we obtain that the 
entanglement entropy scales as \og£ for large £, and the 
effective central charge (calculated for two symmetrically 
placed defects) is given in Eq. (fT3| . 
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